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Quantum tunneling across double potential barriers is studied. With the assumption 
that the real space is a continuum, it is rigorously proved that large barriers of 
arbitrary shapes can be penetrated by low-energy particles with a probability of unity, 
i.e., realization of resonant tunneling (RT), by simply tuning the inter-barrier spacing. 
The results are demonstrated by tunneling of electrons and protons, in which resonant 
and sequential tunneling are distinguished. The critical dependence of tunneling 
probabilities on the barrier positions not only demonstrates the crucial role of phase 
factors, but also points to the possibility of ultrahigh accuracy measurements near 
resonance. By contrast, the existence of a nonzero minimum length puts upper bounds 
on the barrier size and particle mass, beyond which effective RT ceases. A scheme is 
suggested for dealing with the practical difficulties arising from the delocalization of 
particle position due to the uncertainty principle. This work opens a possible avenue 


for experimental tests of the existence of a minimum length based on atomic systems. 


*E-mail: yyanglab @issp.ac.cn 


I. INTRODUCTION 
The scenario of a minimum length (Lmin) plays an essential role in the quantum 
theory of gravity [1-11]. Breakdown of the Lorentz invariance may happen when the 


real space approaches such a minimum length scale [12-20], which is generally taken 


to be the Planck length (lp = ft ~1.6 x 10735 m, fi is the reduced Planck’s constant, 


G is the gravitational constant and c is the speed of light) [7]. Despite numerous 
efforts [1-20], the question remains open regarding the existence of such a minimum 
length [7, 21-23]. Since such a length scale is well below the lower bound of spatial 
resolution achieved by state-of-the-art instruments such as LIGO (~1071? m) [24, 25], 
it is a great challenge for experimental verification. Here, we show the possibility of 
tackling this problem by investigations of quantum tunneling across double potential 
barriers. 

Quantum tunneling [26] is a classically forbidden phenomenon in which a particle 
passes through a potential barrier higher than the energy it possesses. In the early 
years of quantum mechanics, the theories based on quantum tunneling explain some 
puzzles of experimental observations like the thermionic and field-induced emission 
of electrons from metal surfaces [26], and the alpha decay of heavy nuclei [26]. In the 
ensuing decades, researches on the quantum tunneling of electrons in condensed 
matter have led to fruitful discoveries [27-33], and enabled important inventions such 
as the scanning tunneling microscope (STM) [34] and tunneling diodes [35-37]. Since 
the pioneering works by Tsu, Esaki and Chang [31-33], double barriers have received 
a lot of attention while studying electron transport in heterostructures [37, 38]. 
Resonant tunneling (RT) typically takes place in double-barrier systems, in which the 
incident electrons may pass through the barriers without being reflected, i.e., with a 
transmission probability of 100%. Such a behavior is due to the coherent interference 
of electron waves which cancel the reflected waves and enhance the transmitted ones, 
analogous to the resonant transmission through a Fabry-Perot etalon in optics. Typical 
inter-barrier spacing of the devices based on RT is several tens of angstroms (A), 


matching the de Broglie wavelengths of electrons. In recent decades, the phenomena 


2 


of RT in mesoscopic and nanoscale structures continue to attract interests of research 
[39-43]. 

Historically, RT of electrons was considered to gain experimental evidence from 
the negative differential resistance (NDR) found in the current-voltage (Z-V) curves 
[32, 35, 37]. Later, alternative mechanism was suggested for NDR, namely, sequential 
tunneling in which the phase memory of electron wave functions is lost due to 
inelastic scattering [37, 44-51]. It was argued that resonant (coherent) tunneling is a 
prerequisite for sequential tunneling [51]. The effects of external electric field, 
inelastic scattering, and the repulsive interactions between electrons on RT were also 
studied [38, 52, 53]. In spite of these efforts, consensus on the underlying physics is 
yet to be reached. There are still large discrepancies between theory and the measured 
I-V curves (e.g., peak-to-valley ratio). The gap originates partly from the fact that, in 
calculations related to experiments the simplest rectangular barriers (or their variants) 
are adopted, which usually differ significantly from the true barriers felt by electrons. 

To resolve the puzzles, exact theoretical description of the conditions for RT 
across double barriers is highly desired. For the simplest rectangular double barriers, 
exact mathematical relation of energy and geometric conditions has been established 
[38, 54, 55]. For the more general and realistic situation where double barriers are of 
arbitrary shapes, aside from the semi-classical approach [38], full quantum level 
description of the RT conditions is still lacked. It is generally accepted that RT takes 
place when the energy of incident particle matches the energy levels of the 
quasi-bound states within the potential well in-between the two barriers [37, 38, 44, 
49, 51, 52]. In principle, this applies to electrons as well as the massive particles like 
protons, atoms and molecules. Recent simulations have shown the RT of H and He 
atoms across small double barriers, with the barrier height Ep ~ 0.2 eV [56, 57] and ~ 
0.02 eV [58], respectively. However, when a particle tunnels across arbitrarily-shaped 
double-barriers, it is unclear how the level-match condition can be reached, and 
rigorous theoretical descriptions remain elusive. 

In this paper, we revisit this topic in double-barrier systems consisting of equal 


barriers of arbitrary geometries. With the assumption of a continuously varied 
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inter-barrier spacing (equivalently, Lmin = 0), it is rigorously proved that quantum 
tunneling through the double-barrier system with a probability of unity can always 
happen (i.e., RT) when the inter-barrier spacing is appropriately chosen. Exact 
mathematical relation for RT is established. At the presence of a nonzero Lmin (Lmin > 
0), the inter-barrier spacing varies discontinuously, which sets upper bounds for the 
barrier heights and particle mass, above which no RT may happen. The results are 
demonstrated by the tunneling of electrons, protons, and some typical bosons. 
Practically possible scheme is therefore provided for experimental tests of the 
existence of a nonzero minimum length. 

The rest of this paper is organized as follows. Section II presents the analytic and 
numerical results on RT, with examples of typical particles like electrons, protons and 
some bosons. The connection between RT and the continuity of real space is revealed. 
The constraints set by the existence of a minimum length, the practical obstacles due 
to the uncertainty principle and plausible solutions are presented. We conclude in Sec. 


II with discussions on the impacts and future opportunities inspired by this work. 


II. RESULTS AND DISCUSSIONS 

We begin in Part A of this section by performing general analysis on the 
transmission properties of quantum particles across double barriers of arbitrary shapes, 
and prove a theorem which establishes the mathematical condition for resonant 
tunneling (RT). Part B provides analyses on two typical models — rectangular and 
parabolic double barriers. The quantum tunneling of electrons and protons are studied 
and compared, with emphasis on the differences between resonant and sequential 
tunneling. Based on the results of Parts A and B, we show in Part C the upper bounds 
of barrier heights for RT set by the Planck length. In Part D, the fundamental limits 
put by the uncertainty principle are studied and possible solution to position 


delocalization of the incident particles is suggested. 
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FIG. 1. Schematics of resonant tunneling (RT) across double barriers. (a) Quantum 
interference of the incident and reflected matter waves; (b) Modulation of the energy 


levels of the quasi-bound states in-between the two barriers, by varying the 


inter-barrier spacing w; (c) RT spectrum as a function of w with a period of A (= 5). 


A. GENERAL ANALYSES ON ARBITARY DOUBLE BARRIERS 
Generally, double-barriers consist of two identical or different single barriers, 
which are respectively referred to as homo-structured and hetero-structured hereafter. 
The double-barrier considered here is homo-structured in one-dimensional space, as 


shown in Fig. 1(a), with a barrier height E, and barrier width a for each. Our analyses 
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are based on the transfer matrix method, a powerful technique for studying the 
transmission properties in finite systems [31, 59-62]. For the propagation of a 
quantum particle across a single barrier V(x), the transmitted and reflected amplitudes 
(Az, Br; Ar, Br) of the wave functions (W,, Wp) may be related by a transfer matrix 
(denoted by M) as follows [38, 56, 59-62]: 
Ar\ _ A) — (m1 Maiz) (a, 
ay) =M (2) ~ a i) i) () 

The incoming wave function (with incident energy E) is Y, = A,e‘** + B,e~"**, and 
the outgoing wave function is pp = Ape + Bre, where k = ./2mE/h?, m is 


the particle mass. The determinant |M|=1, for systems where time-reversal 


. . i 25 8 be H a 
symmetry preserves, and the transmission coefficient is given by [56] T = mal? = 
11 


mole" In general, the matrix elements m;; (i, j = 1, 2) are complex numbers and 
22 


obey the conjugate relations [56, 59-62] of mj, = M22, and mj, = m31. For a 
homo-structured double-barrier with an inter-barrier spacing w, the following theorem 
holds: 


Theorem. — For any E < Ep», the transmission coefficient (tunneling probability) 


F nt m+0+2ka 
across a homo-structured double-barrier Tpg(E;w) = 1 at w = wp = ae 


where 0 = arg(m,,7), n (referred to as resonance number) belongs to integers. 
Proof. — The updated transfer matrix for a single barrier V(x) translated by a 
distance L = a + w, V(x-L), is given by [57, 63] 
re ma my e7i2kE _ Mas my e72k(a+w) -D 
M21e! M22 migja er M22 


The transfer matrix for the double-barrier (U(x) = V(x) + V(x-L)) is therefore [56] 


(3) 


Mpg = M(L) * M = ( m1 | ie ies 


i2k(a+w m m 
m,e Otw) M22 21 22 


Dn 


The diagonal matrix element describing the transmission properties, (Mpg)11, i 
explicitly calculated to be (Mpg)11 = M112 + mame PKO), Let Z = m? = 
|Zle", ġ = 2k(a +w), mi2m31 = |m,2|* = R, with i the imaginary unit, and the 
angle 0 =arg(Z); the determinant, |M| = 1 = |m,,|* —|m,2|? =|Z| —R_ gives 
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that |Z} =1+R, then (Mpg) = (1+ Roe” + Re~@ = et? + Rel (e-HOt8) 4 
1). When e~'+® = —1, i.e, 6+ 6 = (2n — 1)r, with n being integers, one has 
1 


(Mpp)11 = e. It follows that the transmission coefficient Tpg (E; w) = ene 
DB/11 


1, which corresponds to RT. Using the condition @ + 0 = (2n—1)z, one has 


2k(a+w)+0=(2n—1)m, and consequently, w = = = a= 


= Wp. This 
completes the proof of the theorem. 

It should be stressed here that the proof inherently includes the precondition that 
the inter-barrier spacing w varies continuously (Lmin = 0) such that the angle @ can 
have any desired values to satisfy the RT condition. The theorem points to the 
possibility of penetration of arbitrarily large (but finite) potential barriers by 
low-energy particles with a probability of unity. For a quantum particle with incident 
energy E, it can completely tunnel across a homo-structured double barrier of height 
Ep when the inter-barrier spacing equals w, described above, even in the case E << Ep. 
In addition, one sees that the barrier-barrier separations (w,) for RT are solely 
determined by the parameters (0, a) describing the transmission of single barriers. 
Physically, the onset of RT is due to the presence of quasi-bound states in-between the 
two barriers whose energy levels match that of the incident particles [37, 38, 44, 49, 
51, 52]. A direct consequence is that, any quasi-bound energy levels (E < Ep) can be 
realized within the potential well set by the two barriers via simply tuning the 
inter-barrier spacing, as illustrated in Fig. 1(b). Moreover, from its mathematical 


expression, one sees that (Mpg)11 is the periodic function of w, with a period of 
T= = For a fixed E, the tunneling probability T[E; w] displays periodic variations 
with w, showing comb-like structures with the resonance peaks positioned at L, = a + 
Wn, and the distance between any two neighboring peaks is A= w, — Wn-1 = 7 (Fig. 


1(c)). The value of A is just half the de Broglie wavelength of the incident matter 
wave, indicating the key role of phase factor and quantum interference. Finally, the 
mathematical expression of w, implies that there could be infinitely many resonance 


peaks in free space. The results may be readily extended to two- or three-dimensional 
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systems in that the interaction potentials along the direction of propagation are 


equivalently described by some effective double barriers. 


B. TUNNELING ACROSS TYPICAL DOUBLE BARRIERS: RESONANT 
TUNNELING VERSUS SEQUENTIAL TUNNELING 
For homo-structured rectangular double barriers, analytic expressions of w, are 
available, which enable in-depth understanding of the physics of RT. The matrix 
element mıı describing the transmission across single rectangular barrier (barrier 
height Vo) may be expressed as follows (Appendix A): 
mı = 2ye ~t [i(k? — B?)sinh(Ba) + 2Bkcosh(Ba)], (4) 


1 


where k = §2mE/h?, B =. 2m(V) — E)/h?, y= ape Eq. (4) may be reduced to 
ay = 2ye x gett = 2yoetl ko, (5) 


where o = VA? + B?, A = (k? — B7)sinh(Pa), B = 2Bkcosh(Ba), and the angle 


_ Ay Therefi 2 _ (A+B?) | i2(a—ka) trs; 
a= arctan). Therefore, mj, = “ape © . Using the theorem stated above, 


the angle 6 = 2(a — ka), and then 0 + 2ka = 2a. The inter-barrier spacing is given 


b We T+2Q 
Y Wa > 2k ` 


It follows that 2kw, = (2n — 1)m — 2a, and one arrives at the 


dtanh(Pa) 
1-25?tanh? (Ba) 


B 


equality: tan(2kw,) = , where ô = (E — 5) . Alternatively, this 


equality can be obtained by direct calculation of the squared norm of diagonal 
element |(Mpg)11|?, a function of inter-barrier spacing w: The minimum of 
|(Mpg)11|? leads to RT (Appendix B). The equality for tan(2kw„) is in line with 
Ref. [55], which was derived in a different way. In the special case when the incident 


energy is half the barrier height (E = 0.5Vo), 6 = k, the angle æ = 0, one obtains a 


simplified relation that 2kw, = (2n — 1)r, and w, = Cam = (n — 2) (2), with 


2 
Àa = = is the de Broglie wavelength. In another special case when k <« f and 


pa > 1, i.e., the incident energy is far below the barrier height, one has a = = + 


k k nT 1 ‘ š 7 
F and kw, = nn -zp Wn = op In both situations, the value of w, is 


independent of the barrier width. 
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FIG. 2. Schematics of rectangular (a) and parabolic double barriers (b). The 
inter-barrier spacing of RT (w,), as a function of resonance number (npr), for electron 


(c) and proton (d) across the double barriers, at incident energy E = 0.5 eV. 


The results are applicable to electrons and other massive quantum particles. 
However, demonstration of the quantum interference effects leading to RT would be 
much more challenging due to the large difference in particle masses and the 
corresponding de Broglie’s wavelengths. Here, we perform systematic investigations 
on the RT characteristics of electrons and protons through two model systems: 


rectangular and parabolic double barriers (Fig. 2). Compared to the analytic 
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expressions for rectangular barriers, the transfer matrices for parabolic barriers are 
evaluated numerically [56, 57]. Figures 2(c)-(d) shows the calculated w, for the RT of 
electrons and protons across the two types of double barriers, as a function of the 
resonance numbers (ngr). For the same ner, the w, of electrons is much larger than 
that of proton owing to smaller mass. The different geometries of the potential 
barriers (rectangular vs parabolic) are reflected by the slight differences of w,. Despite 
the differences, the overall comparable magnitudes of the two sets of w, indicate that 
rectangular double barriers may serve as approximations for qualitative description of 
some smoothly varying double barriers with regular geometries. 

At fixed energy E, the tunneling probability varies periodically with inter-barrier 
spacing w. We have further studied such characteristics in case of electrons and 
protons tunneling through rectangular double barriers. Figures 3(a-b) show the 
transmission of electrons, at varying w for E = 0.03 eV and 0.5 eV. The effects of 


incident energy on the tunneling spectrum, T(E; w), are clearly seen. Higher energy 
not only results in smaller period of oscillation (Tt = = but also smaller 


peak-to-valley ratio. The resonance number can extend to very large integers, as long 
as the perturbation from the environment is negligible and the coherence of wave 
functions is maintained. To show the role of coherence, we have studied the 
energy-dependent tunneling probability P(E) of electrons at a fixed inter-barrier 
spacing (w ~ 10um). For resonant (coherent) tunneling, the quantity P(E) (= T(E; w)) 
drops quickly with small deviations from the resonant energy level, Err. For 
sequential tunneling, in which the phase coherence is destroyed in a two-step 
tunneling process, the quantity P(E) is simply product of the transmission coefficient 
across each single barriers: P(E) = T,(E) x T,(E) = 7?(E). Around the resonant 
energy Err (Figs. 3(c-d)), the probability of sequential tunneling (Psr) changes 


smoothly with energy, and is significantly smaller than unity at low incident energies. 
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FIG. 3. Tunneling spectrum of electrons across the rectangular double-barrier shown 
in Fig. 2(a). RT at resonance level Err = 0.03 eV (a) and 0.5 eV (b), as a function of 
inter-barrier spacing w. Panels (c-d): Energy dependence of tunneling probability at 
fixed w, around Err = 0.03 eV (panel c, w = 100008.49 A) and 0.5 eV (panel d, w = 
100991.45 A). The data lines labeled by Pı, Psr and Prr correspond to tunneling 
through single barrier, sequential and resonant tunneling through double-barrier, 


respectively. 


For protons, more radical differences encounter. Shown in Fig. 4(a), is the 
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tunneling spectrum of protons at E = E,/2 = 0.5 eV. The periodically repeated isolated 

lines imply much narrower resonant peaks with comparison to electrons. The enlarged 

structures of one resonant peak are shown in Fig. 4(b). Around the RT peaks, the 

squared norm of transfer matrix element may be expressed as follows (Appendix C): 
|(Mpg)11|2 = 1 + sinh? (2ka) x (kAw)? = 1 + A|M,,|%, . (6) 

where Aw is the deviation from the peak position w,. When AļM11lĝw = 1, T(E; w) 


= 0.5, and one has 


1 


[Aw] = ksinh(2ka) ` 


(7) 

It follows that the term 2|Aw| is the full width at half maximum (FWHM) of 
the resonant peaks. For the double-barrier (Fig. 2(a)) considered here, it turns out that 
|Aw| = 4.235 x 10~15A. When the deviation Aw ~ 10°? A, the tunneling probability 
drops quickly to T(E; w) ~ 10°, in good agreement with the results presented in Fig. 
4(b). In general, given that w is an approximate value to w,, one can determine the 
significant digits of w by designating a deviation Aw when |w — w,| < Aw such that 
T(E;w) > 1 — 8P, where 5P (0 < 6P < 1) is the tolerance of decrease in tunneling 
probability at which significant tunneling (effective RT, events measurable in 
experiment) is maintained. Furthermore, using the proof of the theorem, we find that 
for arbitrarily homo-structured double barriers, the deviation Aw at the tolerance &P 


may be given by (Appendix D): 


i | a SP 
Aw = RG4R) IOP ° (8) 


where R = |m,,|?, k and 6P are defined as above. Here we focus on the case that 
65P = 0.5, which yields the FWHM ( = 2Aw). 

Such ultrahigh sensitivity on tunneling parameters is also found for the RT 
energies. Figure 4(c) compares the tunneling of protons across single and double 
barriers at a fixed inter-barrier spacing (w ~ 20 A). Near resonance, P(E) descends 
drastically from 1 to ~10°'' by a tiny shift of £ = 10" eV from Epr. At the vicinity of 
Err, the dependence of |(Mpg)11|? with deviation AE is given by (Appendix C): 

AE 


2 2 
|(Mpg)1112 = 1 + sinh?(2ka) x (=) x (=) =1+A|M,,|2,. The FWHM at 


12 


energy scale is therefore obtained when A|M,,|i- = 1, and |= =| = ar In 


our case, ~ 4.235 x 10716. When the energy broadening AE = £ = 107 1°eV, 


FE] = 2x 10-2, |(Mpa)ail? = 2.23 x 101, P(E) = |(Mpp)ail"? = 10-135, 


compares well with the numerical results. Without resonance, the probability of a 
two-step tunneling, i.e., sequential tunneling decreases by more than 25 orders of 
magnitude (Fig. 4(c)). The sharp contrast distinguishes RT from sequential tunneling. 
For the more generalized case, the allowed energy broadening AF may be calculated 


as follows (Appendix D): 


(9) 


ee 1 ôP 
= Ga R(1+R) ` 1-6P 


In the case of R >> 1 (large reflection), for instance, tunneling through large barriers 


: ; : ; AE] 1 EA ôP 
or tunneling by massive particles, Eq. (9) is reduced to =| = Pawn nae For a 


single barrier V(x), the reflection and tunneling probabilities (Appendix D) are related 
to R by |r|? = R|t|? = — subjected to the condition |r|? + |t|? = 1. It is 
i 


straightforward that R = — — 1 = — = —_ when R >> 1, where T} (E) is the 
F ltl? T(E) 


tunneling probability across a single barrier. 

As seen from Fig. 4(c), at the absence of phase coherence, the incident protons 
will be nearly completely reflected by a single barrier. On the contrary, when the 
inter-barrier spacing equals w, and phase coherence is maintained, the protons 
penetrate the two barriers with a probability of unity. Such effects are schematically 
illustrated in Fig. 5. The key role of quantum interference is demonstrated. 
Experimental verifications may be carried out using atomically thin membranes, 
which have potential applications as proton sieve filters. Generally, the variation step 
(Al) of the inter-barrier spacing w, required by RT should be the order of magnitude of 
Aw studied above, and no less than the minimum length (1.e., Al ~ Aw > Lmin) such 
that effective RT can be reached by tuning the inter-barrier spacing. This is the topic 


of next subsection. 
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FIG. 4. Tunneling spectrum of protons across the rectangular double-barrier shown in 
Fig.2(a). Resonance at Err = 0.5 eV (a). Panels (b-c): Variations of tunneling 
probability with respect to small deviations from the RT parameters: (b) Inter-barrier 
spacing w = (Ny - Np)XAl + Wn, Ny = 596 and Al = 10°'° A; (c) Incident energy at the 
vicinity of Err, for w = 20.137016632763302 A, and the energy deviation ¢ = 10" eV. 


All digits of w are meaningful. 
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FIG. 5. Schematics of quantum tunneling of protons across single barrier (upper panel) 
and double barriers (lower panel) at the presence of RT. The probability of reflection 


is denoted by Pr; and tunneling by Pr, i = 1, 2. 


C. UPPER BOUNDS OF RT BARRIERS SET BY THE PLANCK LENGTH 

The critical dependence of the tunneling probabilities on the barrier positions not 
only demonstrates the crucial role of phase factors, but also points to the possibility of 
ultrahigh accuracy measurements near resonance. As shown above, the deviation of 
|Aw| = 4.235 x 10715A leads to a 50% drop of P(E) of protons across a rectangular 
double-barrier. Such a deviation is several orders of magnitude below the smallest 
length scale sensed by LIGO [24, 25]. Even smaller |Aw| is expected for heavier 
particles or larger barriers. As mentioned above, to have measurable RT within some 
tolerance 6P, an upper bound of deviation from the exact peak position w, is given by 
Aw = |w — w,|. Suppose that the elementary variation step of distance is Al; if the 
real space is a continuum (Lmin = 0), then Al can be arbitrarily small and in principle 
w, can always be reached by a finite number of operations. In this case, the theorem 
stated above always holds. On the contrary, if a nonzero Lmin exits (Lmin > 0), to have 
significant RT, the variation step should satisfy Aw > Al > Lmin. Let n = [Log 9(Aw)], 
then it is a feasible choice to set Al = 10”, such that w = NxAI, with N being an 


integer and the modulo w mod Al = 0. Consequently, the existence of a minimum 
15 


length leads to the inequality of |Aw| > Lmin, and therefore puts some upper bounds 
for the particle mass, barrier height and barrier width, above which RT will cease. In 


the special case of tunneling across rectangular double barriers at E = 0.5E,, 


i : . 
> Lmin, Which may be rewritten as 


realization of RT requires that |Aw| = ksinh(2ka) — 


xsinh(y) <—, (10) 


Lmin 


where x = 2ka, k = ee = =. The upper bound of the term mVp is 


therefore determined for a given barrier with a. Assuming that the minimum length is 


identical to the Planck length (Lmin = lp), the upper bounds of barrier height (Vmax) 
for electrons and protons are calculated and shown in Fig. 6. At a barrier width of a = 
1 A, 5 A, and 10 A, the Vmax is ~ 5652.72 eV, 239.42 eV, and 61.32 eV for electrons, 
and is ~ 3.08 eV, 0.13 eV, and 0.03 eV for protons, respectively. It is seen that Vinax of 
RT decreases fast with increasing barrier width. For instance, when the barrier width 
increases to a = 20 and 30 A, the value of Vmax is respectively ~ 15.70 and 7.08 eV for 
electrons, which may be feasible for experimental tests using metal-insulator-metal 
double barriers. Meanwhile, the same variation trend is found for electrons and 
Mme _ 1 


protons, with the magnitude of Vmax being scaled by a factor of ņn = a Ge where 
p 


Me and mp is respectively the mass of electron and proton. This is due to the conjugate 
relation that the particle mass times barrier height (mV ) is a constant at fixed barrier 
width. For the general case of tunneling through arbitrary double barriers, the 
constraint imposed on the particle mass and barrier size due to a nonzero minimum 
length is given by (Appendix D): 


h 1 ôP 


2V2mE «| R(1+R) ap = Lmin > (11) 


where R = |m,,|?, and &P (0 < 6P < 1) has the same meaning as above. Provided 


that Lmin = lp and ôP = 0.5, the inequality reduces to 


i 1 
2V2mE 4| R(1+R) 2 lp . (12) 
Since the parameter R is generally an increasing function of barrier size (Appendix D), 
the upper bounds on the barrier size of RT are therefore determined by the Planck 


length. 
16 


electron 


0.20 


0.15 


0.10 


(keV) 


0.05 


max 


V 


0.00 
6 


0 2 4 6 8 10 
Barrier Width (A) 


proton 


0.15 
0.12 
0.09 
0.06 


(eV) 


0.03 
0.00 


max 


V 


0 2 4 6 8 10 
Barrier Width (A) 


FIG. 6. Calculated upper bounds (Vmax) of the barrier height of rectangular double 
barriers set by the Planck length for electrons (a) and protons (b), as a function of 


barrier width. The values of Vmax at barrier width of 6-10 A are highlighted in the 


insets. 


D. FUNDAMENTAL LIMITS PUT BY THE UNCERTAINTY PRINCIPLE 
AND POSSIBLE SOLUTION 


For a group of incident particles, given that the standard derivation of energy 
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distribution, og, is approximately the term AE for A|M,,|%; = 1 and P(E) = 0.5. The 
narrow window of energy dispersion implies that the momenta of the particles 
distribute dominantly within a narrow interval with a small standard derivation (Ap). 
As a consequence of the uncertainty principle, the standard derivation of position, Ax, 
is expected to be large. Table I lists the energy and momentum broadening, and 
estimated standard position derivations of protons when P(E) = 0.5 at E = 0.5E;, fora 
number of rectangular double barriers with Æ, = 1, 0.5, 0.2, and 0.1 eV, and w ~ 20 A. 
It is clearly seen that the energy broadening increases significantly with decreasing 
barrier height, resulting in reduced standard derivations of position. For Ep = 1 eV, the 
requirement of ultrahigh monochromaticity of incident energies leads to a very small 
Ap and consequently a quite large Ax (~ 1.53x10* m), which is practically very 
challenging, if not impossible for experimental tests. Much smaller Ax (9.25x10° m) 
is found when E> decreases to 0.1 eV. 

It should be stressed here that, the constraint on the standard deviations of 
particle momentum and position imposed by the uncertainty principle does not 
exclude the possibility that a subgroup of particles with ultrahigh monochromaticity 
coexist with another subgroup of particles with large energy broadening. The reason is 
that, by definition, the standard derivation of some physical quantity of a single 
particle is the statistical average over a large number of events, which may be 
equivalently evaluated by the statistical results of a large number of identical particles 
within a small interval of time. Indeed, this is in line with the impossibility of 
measuring the quantum state of a single system [64]. Therefore, a possible recipe to 
the practical difficulty of position delocalization is to have much larger standard 
derivation of kinetic energy distribution than that required for the half drop of P(E), 
i.e., og >>AE, such that the standard derivation of momentum, o,, is much larger than 
the Ap corresponding to AE. Consider two microcanonical ensembles containing N 


weakly interacting identical bosons which follow the kinetic energy distributions g)(E) 


and g:(E) respectively: N = if, gı (E)dE = J, g2 (E)dE. In addition, they have the 


same averaged kinetic energies: (E) = i Eg, (E)dE = i Eg, (E)dE. The key 
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difference is the standard deviation of kinetic energies: Og > Ogpı~AE, i.e., the 
energy broadening of the first group of particles (distribution described by g;(E)) is 
approximately the energy deviation for the half drop of P(E), while is much smaller 


than that of the second group. The distribution function of the mixed 2N-particle 


2 2 
ensemble is g(E) = g,(F) + g2(E), with the standard deviation og = [otot z 


z >> AE. Therefore, mixing of the two groups of identical particles has drastically 


increased the energy broadening and reduced the position uncertainty. Meanwhile, 
sufficient number of particles for resonant transmission is maintained. The 
modifications introduced by the procedure are illustrated in Fig. 7. For the general 
case of P(E) = 1 — òP, the energy broadening AE is given by Eq. (9) and can be 
similarly analyzed. In weakly interacting dilute atomic gases, two-body collisions 
dominate the interactions which simply exchange particle momenta and therefore 
keep the kinetic energy distributions unchanged. 

In practice, the first group of particles may be prepared using the Bose-Einstein 
condensates [65-67], in which the momenta of all involved bosons are expected to 
have approximately the same value: Condensation in the momentum space. The 
second group of particles may be prepared at temperatures slightly above the critical 
temperature T. of phase transition from normal states to the new quantum states like 
superconductivity, superfluidity, or Bose-Einstein condensation. As an example, the 
RT of some typical bosons (Cooper pair of superconducting Nb, “He, "Li, Na, *’Rb) 
across rectangular double barriers is studied and the related parameters are presented 
in Table II. The effects of energy broadening through mixing identical bosons of 
different ensembles are evidenced by the significantly reduced standard position 
deviations. Nevertheless, preparation of the first group of particles with ultrahigh 
monochromaticity remains challenging even with state-of-the-art technique. Another 
challenge to experimental tests may be the acceleration of the condensates as a whole 


to desired incident velocities while maintaining the states of condensation [68, 69]. 
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AE) 


Energy 


FIG. 7. Schematic diagram for the kinetic energy distribution (f{(E)) of identical 
particles in microcanonical ensembles: Particle groups of high monochromaticity (E) 
= g1(E)), low monochromaticity (f(E) = go(E)), and their superposition (KE) = g(E) = 
gi(E) + g2(E)). 


Table I. Parameters describing the RT of protons across rectangular double barriers at 
E = 0.5E,. With the deviation of AE or |Aw| from the parameters for resonance, the 
tunneling probability drops from 1 to 0.5. The corresponding momentum broadening 
Ap, and the minimum standard deviation of particle positions Ax,, are calculated using 


the relation ApAX» > f/2. In all cases the barrier width a = 1 A. 


E, (eV) w (A) AE (eV) |Aw|(A) | Ap (kg.m/s) AXm (m) 
1 20.137016632763302 | 2.103x10'° | 4.235x107° | 3.443x10°° |  15314.7 
0.5 20.17790917547 1.320x10 1 | 5.328x10"' | 3.056x10~° 1.725 
0.2 20.1380336 2.671x10° | 2.690x107 | 9.778x10° | 5.392x10 
0.1 20.15963 1.101x107 | 2.220x10° | 5.700x10°° | 9.251x10° 
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Table II. Similar to Table I but for the RT of some typical bosons with incident 
energy E = 0.5E,. In all cases the barrier width a = 1 A, and barrier height E, = 
0.01 Vmax, With Vmax being the upper bound set by the Planck length. The Cooper pairs 
of electrons are represented by e...e. The energy broadening and resulted 


uncertainties of momenta and positions of mixed particle groups are displayed in the 


lower lines of the same columns. The broadening parameter of energy is chosen such 


that og = kpT., with kg the Boltzmann constant and T, the phase transition 


temperatures. 
: , AE (eV Ap (kg.m/ AXm 
Boson | £E, (eV) w (A) \Aw| (A) NE APR i 
og(eV) 6, (kg.m/s) 6, (m) 
e.e , | 141x10'? | 1.43x10°7 368.15 
28.26 6.3439 3.16x10™ ; = - 
(in Nb) 1x10 1.02x10 5.19x107 
P 5 a | 3.83x107° | 1.43x10°7 368.15 
He 7.69x10 6.3439 3.16x10 : 5 > 
5x10 1.87x10 | 2.82x10" 
n ‘ a | 219x10" | 1.43x10°7 368.15 
Li 4.39x10 6.3439 3.16x10 
1x10? 6.54x10° | 8.07x10* 
5 P a | 6.67x10" | 1.43x10°’ 368.15 
Na 1.34x10 6.3439 3.16x10 
1x10? 2.15x10°! | 2.45x104 
a à a | 176x10" | 1.43x10°’ 368.15 
Rb | 3.53x10 6.3439 3.16x10” 
1x10? 8.13x10°! | 6.49x10° 
HI. CONCLUSIONS 


To summarize, we have studied quantum tunneling across double barriers and 
arrived at a theorem which leads to several physical consequences. First of all, by 
tuning the inter-barrier spacing, it is possible that low-energy particles penetrate 
arbitrary finite-sized potential barriers completely via resonant tunneling (RT). This 
result points to the possibility of significant tunneling of massive quantum particles 
across large barriers at mild conditions. Secondly, it is possible to construct any 
desired quasi-bound energy levels within the quantum well formed by the two barriers 
via adjustment of the inter-barrier spacing. Thirdly, for the RT of quantum particles, it 


is possible to detect the tiny variations of energy levels and positions of the involved 
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potential barriers with unprecedented accuracies. Finally, the critical dependence on 
inter-barrier spacing (consequently the phase difference) demonstrates again the vital 
role of phase factor of wave function, which has manifested itself in some remarkable 
phenomenon such as the Aharonov-Bohm effect [70]. 

Demonstration of the above mentioned results involves two key factors: (i) 
Continuity of the real space and (ii) Energy monochromaticity of the incident particles. 
The first is determined by whether or not a nonzero minimum length (Lmin) exists, and 
the second is affected by the uncertainty principle. Provided that Lmin = 0, the 
distances in real space change continuously and RT can always be realized at given 
incident energies. On the contrary, the existence of a nonzero Lmin will set constraints 
(upper bounds) for the particle mass, barrier height, and barrier width, beyond which 
no RT is expected. Meanwhile, to surmount the practical difficulty (position 
delocalization of incident particles) owing to the uncertainty principle, we suggest a 
plausible scheme in which the high- and low-monochromatic flows of identical 
particle groups are mixed. Potential applications of Bose-Einstein condensates in the 
scheme are discussed. This work reveals the deep connection between two seemingly 
different branches of quantum physics: quantum tunneling and quantum gravity, and 


opens a possible avenue for testing the existence of a minimum length. 
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APPENDIX A: MATRIX ELEMENT my FOR TUNNELING ACROSS SINGLE 
RECTANGULAR BARRIER 

Within the transfer matrix method, we derive the diagonal matrix element mıı that 

describes the transmission across a single rectangular barrier. For a quantum particle 

with incident energy E tunneling through a rectangular barrier with the height of Vo 


and width a, the transfer matrix may be given by [56]: 
y = 2 (Gk+ Be Pe (ie — poe ere a +ik B- a 
1T Zik\ (ik — Beh (ik + BeBe JZB\B -ik B+ ik 
y (ik + B)2e~Gk-B)a — (ik — Bye Get Aa (B? + k?)(e~tk-Ba = e~(ik+B)a) 
== (p2 + k?)(e Pe _ elik+b)a) (ik + p)2elik-P)a — (ik — Bye kt Ba 


1 
where k = /2mE/h?, B = [2m(Vo— E)/Ñ?, y = zi 


The first diagonal term is: m4, = ” [Cik + B)2e-Gk-B)a — (ik — p)%e~Gk+B)2] = 
~iy[(ik + B)2e~U- 4 — (ik — B)2e~“**+ 2), which can be reduced to 
—2iye‘**[(B? — k?) sinh(Ba) + 2iBkcosh(Ba)], and finally one has 


2ye “*4\i(k? — £?) sinh(Ba) + 2Bkcosh(Ba)] (Al) 


M11 


M11 


APPENDIX B: DEDUCTION OF ALTERNATIVE RT CONDITION 
In this appendix, we deduce the resonant tunneling (RT) condition for 
homo-structured rectangular double-barriers. 
For a double-barrier (DB) consisting of single rectangular barriers with the height Vo 
and barrier widths a and b, the diagonal element Mj, of the transfer matrix M may be 


expressed as follows [57]: 


|M111? = 
1+ (Bzr) [sinh? (£b) + sinh?(Ba)] + 2 e sinh? (Fb)sinh? ($a) 
4ß?k?2 4ß2k2 
-r + k?)? sinh(Bb) sinh(Ba) { [ [(B? + k2)? — 8k? ß?] cosh £ (a + b) 
— (B? + k*)? cosh £ (a — b) ] : cos2[k(b + w) — ka] 
—4kB(B? — k?) sinh f (a + b) : sin2[k(b + w) — ka] }. (B1) 


where k = y42mE/Ř?, B = 42m(Vo — E)/h?, E is energy of the incident particle. 
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In the case of homo-structured DB, a = b, then 


Mal? = 1+ FED x [asinn?(pay] +2 [E] sinh* a) - (6? + 
k?)? sinh?(Ba) x {[ [(B2 + k”)? — 8k? ß?] cosh(2Ba) — (B2 + k?)?] cos(2kw) — 
AkB(B? — k?) sinh(2 Ba) - sin(2kw)}. (B2) 


For a given E, |M,,|* is the function of inter-barrier spacing w, the 


= 2 
T(E; w) 
minimum of |M,,|? gives the maximum of transmission coefficient T(E; w), i.e., 
resonant tunneling (RT). The condition of RT can be established by Z [Mial? = 0. It 
follows that, 
[ [(82 + k2)? — 8k? ß?] cosh(2Ba) — (B2 + k?)?] x (—2k) sin(2kw) — 
AkB(B* — k?) sinh(2 Ba) x (2k) cos(2kw) = 0, and consequently 


2 4kB(B2-k?) sinh(2 Ba) 
tan(2kw) ~ (B2+k2)2-[(B2+k2)2-8k2 B2] cosh(2Ba) (B3) 


By dividing the term $k? in both numerator and denominator, Eq. (B3) changes to 


enChne 4( fp) sinh Ba) _ 46 sinh(2 Ba) _ 
(£ =) -E +4) e] cosh(2pa)  (52+4)-(6?—-4) cosh(2Ba) 
ô sinh(2 Ba) 


= (E_E 
(1+562)+(1-252) cosh(2Ba)’ WEIS - ) 


Recalling that sinh(2 Ba) = 2 sinh( Ba) cosh( Ba), cosh(2Ba) = 2 cosh?( Ba) — 1, 


one has 
= 26 sinh( Ba) cosh(Ba) __ 26 sinh( Ba) cosh(Ba) i 
muh) = (1+762)+(1-282)(2cosh2(Ba)-1) — 35?+(1-6?)xcosh?( fa) ” wee pai UE 
reduced to 
tan(2kw) = U _ _S tant) (B4) 
(a) T (=h z sech? (Ba) 


Using the equality sech?(Ba) = 1 — tanh? (£a), Eq.(B4) is finally reduced to 
6 tanh( Pa) 


62 
1—-—-tanh? (Ba) 


tan(2kw) = (BS) 


APPENDIX C: DEPENDENCE OF TUNNELING ON SMALL POSITION AND 
ENERGY CHANGES 
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In this appendix, we deduce the mathematical expressions describing the 
dependence of squared norm of diagonal transfer matrix element, |(Mpz)41|7, with 
respect to slight deviations from the peak positions and incident energies at resonant 
tunneling (RT), for the special case when the incident energy is half the barrier height 
(Vo) of a homo-structured rectangular double-barrier (width of single barrier: a). The 
inverse of |(Mpg)11|2 then describes the dependence of tunneling behavior on small 
position and energy changes. 

In general, |(Mpg)11|* = f(E;w), is the function of incident energy E and 
inter-barrier spacing w. At the vicinity of RT, the function f (E;w) can be expressed 
as functions of small deviations from RT parameters using the Taylor series, by 


considering the fact that |(Mpg)11|? = 1 and (=) = 0, (2) = 0 at the RT point. 


I. For constant E, the dependence on deviation (Aw) from the RT positions (w,) is 
1 (a? 
|(Mpp)iil? = fŒ; w) = 1+5(54) x (aw)? 214+ AlMuli, = (CD 


Using the expressions for rectangular double barriers (Appendix B), one has 


of = (+) (6? + k?)*sinh?(Ba)[g(B, k) x (2k) x sin(2kw) + h(B, K) x 
(2k) x cos(2kw)], (C2) 


af _ (B?+k?) sinh? (pa) 
ðw2 AB*k2 


where k = /2mE/h?, P =./2m(V) — E)/h?, 
g(B,k) = [(B? + k”)? — 8B7k*] cosh(2fa) — (B? + k*)* , h(B,k) = —4B (8? — 
k?)sinh(2Ba). 


[g(G, k) cos(2kw) + h(B, k)sin(2kw)], (C3) 


The condition (5) = 0 gives that g(6,k) sin(2kw) = h(8,k)cos(2kw), and then 


_ h(B,k) 
tan(2kw) = ae (C4) 
2 2442) oi nn2 
Using Eq. (C4), <4 = en heg, k)sin(2kw){cot2(2kw) + 1], and then 
a2f _ (B2+k?)’sinh2(Ba) — h(B,k) 
aw2 4B4k2 Sati) (C5) 


For rectangular double barriers, we have the general relation 2kw = (2n — 1)r — 2a, 


2_ a2 
and @ = arctan (KS ) tanh(£a)]. 
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Consequently, 


(k?-6?) 
2tana a, OBR tanh(Pa) 


sin(2kw) = sin(2a@) = 


1+tan2a 1+tan2a 


Finally, 


af _ (B2+k?)’sinh?(Ba) h(Bk) _ (B?+k?) sinh?(2ga) 
ðw? AB*k2 sin(2kw) 22 


(1 + tan?) (C6) 


2 
It is clear that ot > 0 holds for all allowed incident energies E, which proves that 


the term |(Mpg)11|? arrives at its minimum and its reciprocal gives the maximum of 
transmission probability, i.e., 1. 


When the incident energy is half the barrier height, 6 = k, we have a = 0, and 
af 2.. 12 
P 2k*sinh* (2ka), and therefore 


2 
I(Mpp)1|? = 1+ (22) x (Aw)? = 1 + sinh? (2ka) x (kAw)? (C7) 


II. For constant w, the dependence on deviation (AE) from the RT energies (Err) is 


1 (a? 
|(Mpp)iil? = fŒ; w) = 1+5(S4) x (QE)? = 1+ AIM le (C8) 
2 2 
Compared to (22), computation of (£) is much more complicated. Alternatively, 


we directly consider the dependence of |(Mpg)11|7 with energy deviation (AE) to the 
second order. For the special situation 6 =k, the mathematical expression of 
|(Mpg)11|? is reduced to 

|(Mpg)11|? = 1 + 2sinh? (ka) + 2sinh*(ka) + sinh*(ka)[cosh(2ka) + 1] x 
cos(2kw) (C9) 

Recalling that 2kw = (2n — 1)r for 6 =k, the term cos(2kw) may be expressed 


by Taylor series around the RT point with respect to Ak to the second order: 
cos(2kw) = —1 + = (2w)? (Ak)? = —1 + 2(wAk)? (C10) 

Substitution of cos(2kw) with Eq. (C10) leads to the following 
|(Mpg)11|? = 1 + sinh?(2ka)(wAk)? (C11) 


Using k =./2mE/h?2, and then Ak = 2E Nm =*x = one finally arrives at 
8 2VE 2 E 
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|(Mpp)11l? = 1 + sinh?(2ka) 2)? (=)? (C12) 


APPENDIX D: GENERALIZED CONSTRAINTS ON BARRIER SIZE DUE 
TO A MINIMUM LENGTH 
In this appendix, we deduce the constraint on the barrier size (barrier height, 
barrier width) for effective resonant tunneling (RT) at the presence of a nonzero 
minimum length (Lmin). As defined above, effective RT implies that giving the 
deviation Aw when |w — w,,| < Aw, the inequality T(E;w) = 1 — 8P holds, where 
P (0 < P <1) is the tolerance of decrease in tunneling probability at which 


significant tunneling is measurable. Based on the proof of the theorem, we have 


1 1 1 


Tpg(E; w) “ena ee a oe n l (D1) 
Then Aw is determined by the equality as follows 
oe =1-— ðP. (D2) 
Equivalently, 
|1 + RED + YI =. (D3) 
It follows that 
|1 + R + R(cos(¢ + 6) — isin(d + 6))|? = — : (D4) 
(1 +R + R cos(ġ + 0))? + R?sin? ($ +0) = — , (D5) 
(1 +R)? + R? + 2(1 + R)R cos(ġ + 0) = — , (D6) 
and then reduces to 
2R(1 + R) [1 + cos( + 0)] = = . (D7) 
Consequently, we arrive at 
cos(# + 0) = -1 + —— x. (D8) 


2R(1+R) ` 1—-ô6P 
For a homo-structured double-barrier system, the two parameters 0 = arg(m,,7), 
and R = |m;2|? are solely determined by a single barrier V(x). The tunable 


parameter is @ = 2k(a + w), via variation of the inter-barrier spacing w by a small 
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magnitude of Aw. At the vicinity of RT, |Aw| << wn. Around @ + 8 = (2n — 1)r, i.e., 


the RT points, expansion of cos(@ + @) using Taylor series to the 2" order, we have 
cos(@ + 8) = -1 += (Ag)? , (D9) 


where Ad = +2kAw. Comparison of Eq. (D8) and (D9) gives that 


2 ga, (D10) 


2kAw = R(1+R) ^ 1-6P 


Finally, we get 


eo (D11) 


X —— > Lmin (D12) 


For a single barrier V(x), the reflection coefficient is given by [59-62] 


2 _ Iml? _ = 2 
Irl = —*, = RT, (E) = Rel’, (D13) 


[m11]? 
where R = |m;2|?, and T(E) = |t|? is the transmission coefficient across V(x) at 
energy E. Conservation of probability current gives that |r|? + |t|? = 1. Qualitatively, 
|r|? increases with barrier width a and barrier height E,, which indicates that R is the 


increasing function of barrier size parameters a and Ep: R = R (a, Ep). Larger barrier 


size results in larger value of R. Substitution of k with ane the inequality (D12) is 


therefore 


A ae: 
2V2mE «| R(1+R)  1-SP 


= Lmin (D14) 


This is the constraint imposed on the particle mass, barrier height, and barrier width 


due to the minimum length. 


In the case 6P = 0.5, FWHM ( = 2Aw) is obtained. Given that Lmin = lp, we have 


A 1 
2V2mE 4| R(1+R) 


=i. (D15) 


For a fixed particle mass m and incident energy E, the inequality (D15) sets upper 


bounds on R and consequently the upper bounds for barrier size of V(x): the barrier 
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width a and barrier height Ep. 
Furthermore, we can derive the constraint on the broadening of incident energy 


by using Eq. (D9). In this case, w (= w,) and a are fixed, Ah = 2Ak(a + w). Using 


k= an we have Ak =ex then 2Ak =kx =, and Ad = k(a +w) x=. It 
follows that 
Fe ul ôP 
(Ag)" = R(1+R) Xip ` (D16) 
Then 
AE 1 8P 
Ag = k(a + w) x |=| = RAFR) I-P (D17) 
Finally we have 
AE 1 1 6P 
| = kaw ROIR) 1-5P (D18) 
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